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1. Introduction 

Consider a feedback control loop with its inputs and its outputs partitioned in a compatible 
way: Mi =(mh,«i 2 , • • • ,Mi„) , u 2 = (u 2l ,u 22 , ■ ' ■ ,«2«) . Ji =(y 11 .^ 12 , •” »yii») . and 
y 2 = (J21 >y 22 » ‘ ' ' iy\n) T ■ The controller C is decentralized iff it is block diagonal i.e. the i-th 
subvector, y of the manipulated variable vector is only affected by the i—th subvector, y^, of 
the measured output: 

y u = Cii(uu-y 2i ) 

Since the early 70’ s significant research effort has been expended on the subject of decen- 
tralized control. Nevertheless, two unanswered questions remain : 

(a) given the set of measurements and manipulations, how does one select the appropriate pair- 
ings ? 

(b) How can one assess fundamental limitations to decentralized control system performance ? 

The first question is referred to as the decentralized control structure synthesis problem 
while the second can be unequivocally addressed only through the optimal decentralized con- 
troller synthesis problem. Given the set of measurements and manipulations, the solution of the 
decentralized controller structure synthesis problem determines the flow of information in the 
control loop, or equivalently the pairings between the measurements and the manipulations. The 
solution to the second problem determines the best achievable closed-loop dynamic performance 
for the given decentralized control structure. 

It has been established, that given a plant and a decentralized structure there may not exist 
any decentralized stabilizing controllers with that structure. Aoki (1972) [2], demonstrated that 
there may exist decentralized control structures that prevent stabilization of the closed loop. 
Wang & Davison (1973) f 16], introduced the notion of decentralized fixed eigenvalues also 
called as fixed modes of a given system. Algebraic characterizations of the notion of decentral- 
ized controllability, which is related to the fixed mode concept, for the two input vector case are 
given in Morse (1973) [1 1], Corfmat & Morse (1976) [3], [4], and Potter, Anderson & Morse 
(1979) [12]. Anderson & Clements (1981) [1], employed algebraic concepts and characterized 
the decentralized fixed eigenvalues of a system and presented computational tests for the 
existence of fixed modes. 

Recently, the issue of stability of decentralized control systems has been addressed within 
the fractional representation approach to control theory. For linear time invariant processes, 
Manousiouthakis (1989) [9] presented a parametrization of all decentralized stabilizing controll- 
ers for a given process and a fixed decentralized control structure. Within this framework, any 
decentralized stabilizing controller is parametrized in terms of a stable transfer function matrix 
that has to satisfy a finite number of quadratic equality constraints. For the same class of 
processes (LTI plants) Desoer and Gundes presented an equivalent parametrization where the 
stable parameter satisfies a unimodularity condition (Desoer & Gundes, 1990, p.122, 165) [7]. 

In this paper, the Manousiouthakis parametrization is employed in mathematically formu- 
lating the optimal controller synthesis problem. The decentralized performance problem is for- 
mulated as an infinite dimensional /' optimization problem. Performing appropriate truncations 
a finite dimensional optimization problem is obtained. Theorems that establish the connection 
between the two problems are presented. It is shown that iterative solution of the finite dimen- 
sional problem creates a sequence of values that converges to the values of the infinite 
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dimensional problem. Based on these conergence results a computational procedure that yields 
e-optimal solutions to /* optimization problem is outlined. Locally optimal solutions to the 
intermediate finite dimensional problems can be obtained through existing nonlinear optimiza- 
tion algorithms (MINOS, GINO etc.). Global solution of the intermediate finite dimensional 
approximations guarantees that the limit of the sequence that is being created corresponds to the 
best performance that can be obtained by the given decentralized structure. Feasibility (or 
infeasibility) of the optimization problem is equivalent to existence (or nonexistence) of decen- 
tralized controllers with the given structure. 

2. Mathematical Preliminaries 

2.1. Fractional Representations of Linear systems 

Let G be the set of all proper, rational transfer functions and M ( G ) be the set of matrices 
with entries that belong to G. Also let 5 be the set that includes only the stable members of G 
and let M (S) be the set of matrices with entries that belong to S. 

In this work, theoretical results related to the notion of doubly coprime fractional representa- 
tions and the parametrization of all stabilizing compensators are used. These results and a com- 
plete exposition of the underlying theory can be found in Vidyasagar (1985; pp. 79, 83, 108, 

1 10) [14]. The notation used in the present work is compatible with the notation in the 
aforementioned reference. 

2.2. Input - Output Linear Operators 

One of the frameworks developed to describe the stability and performance of dynamical 
systems is the input - output approach. Although the theory has been developed for both continu- 
ous and discrete systems, in this work the focus is on discrete systems. 

In the sequel the fact that every linear BIBO operator can be represented by an l x sequence 
will be utilized. For such operators the /” - /°° induced norm is equal to the / 1 norm of the 
corresponding / 1 sequence. The results that are used can be found in Desoer and Vidyasagar 
(1975; pp. 23-24, 100, 239) [5]. 

2.3. Elements from Real and Functional Analysis 

The notion of denseness will be used in the proofs in Section 4. The fact that the set <j) 0 of 
all sequences with finitely many nonzero elements is dense in /' will also be used. Properties of 
the compact sets will be used in Lemma 2 in Section 4. The related theory is given by Wheeden 
& Zygmud (1977, pp. 4, 8-9, 134) [17]. 

The properties of point-to-set mappings are also used. All relevant results can be found in 
Fiacco (1983; pp. 12, 14) [6]. 

3. Parametrization of Decentralized Stabilizing Controllers 

In this section, the results presented by Manousiouthakis (1989) [9] are outlined. The 2- 
channel case is outlined in section 3.1. In section 3.2, the result corresponing to the general l- 
channel case is presented. 
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3.1. 2-Channel Decentralized Control 

Consider a feedback control system shown with plant P and controller C: 


C = 


Cl,, C ,,i 

C 2, 1 C2,1 


e G nxm , C u e G n ' xm ' , C 2 , 2 e G (n n,)x(m m,) 


P = 


p 1.1 

p 2,\ p 2,\ 


e G mxn , P],i e C m,xn ' ,P 2 ,2 6 C 


(m -m t )x(n -rc j ) 


Manousiouthakis (1989) [9], demonstrated that based on the YJB parametrization of all stabiliz- 
ing compensators for a given plant, the set of all 2-channel decentralized stabilizing compensa- 
tors for the given plant can be described as: 

S d (/>)={ C = (X +DpQ ) (Y -NpQ ) -1 ,det(Y-N P Q)*0,Qe A/(S); 


Si + QS 2 +S 3 Q + QS 4 Q-O, 


( 1 ) 


where, 


5, =YL ni X-XL m Y S 3 =YL ni D P +XL mi N P 
S 2 =N P L n X+D P L mi Y S 4 =NpL ni Dp-DpL m N P 


(la) 


and 




/„ 


0 

0 


= L 


-1 


n , ’ 


L m , — 


f ™> 0 
0 ^ m -m 1 


= L 


-1 


3.2. 1-Channel Decentralized Control 

The parametrization of all Ixl block diagonal stabilizing controllers is based on the results 
of the previous section, namely relations (1), and (la). It has been established that the set of all 
/-channel decentralized stabilizing controllers can be parametrized as (Manousiouthakis, 1989) 
[9]: 




C = (X + DpQ)(Y-N P Qr l , dei(Y-N P Q)*Q , Qe M(5); 


S l j + QS 2 j + S v Q+QS 4 j Q=0J=\, 


( 4 ) 


The transfer matrices S„ , / = 1, 2,3,4 , j = 1, ■ ■ • ,/— 1 are given by relations similar to (la). 


398 



4. Optimal Decentralized Control System Performance 
4.1. I 1 Performance Results 

The performance problem is often posed as follows : determine whether the output of the 
system remains within specified bounds for all bounded external disturbances and for all times. 
The idea of considering disturbances bounded in magnituded was introduced by Vidyasagar 
(1986) [15] and led to the ^-optimal control problem. 

The mathematical formulation of the problem is performed as follows: 

(a) Disturbance, d, is bounded 

0<|d(£)|< w\ , Vk>0<=> |K^|U<i. 

(b) Output, y, is bounded 

0<|y(*)|< W 2 1 , V & > 0 <=> ||»v 2 y lU ^ 1 • 

(c) Satisfy bounds on y for all allowable d 

<D(P,C) = sup ||w 2 y||oo<l. 

||w'7 1 rffloo ^ J 

Let H(P,C) be the closed loop map between the disturbance ( d ) and the output (y). Then, 
y = H(P,C)d. According to (c) the output of the system remains within the desirable bounds iff: 

0(P,C)<1 <=> sup |w 2 y || « = sup ||w 2 //(P,C)w 1 (w7 1 d)||„ = 

IK’rfHooSi IK'4||oo£i 

= || w 2 H(P,C)w l H^^l 

To determine what is the best performance that a decentralized controller can deliver, the 
value of the following optimization problem should be identified : 

v = inf ||w 2 // (P,C)w j ||oooo 
CeSAP ) 

If the value of this problem is less than 1 then there exist controllers with the given structure 
such that the output of the system satisfies the performance requirements. For simplicity in nota- 
tion, w x and w 2 will be augmented in the map H (P,C). 

Employing the parametrization of all stabilizing decentralized controllers the last optimiza- 
tion problem is expressed in terms of the stable map Q e M (S) : 

v = inf |//(P,Q)|| inf ||r, - T 2 Q T 3 1|_ (DPP) 

QeM(S) QeM(S) 

subject to, 

5, +QS 2 + S 2 Q + QS 4 Q=0 
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In the last problem, the optimization variable belongs to an infinite dimensional linear 
space. In addition, the closed loop map H{P,Q) is affine \nQ:H (P,Q)=T j -T 2 QT 3 where 
T j ,T 2 ,T 3 e M (5), are known and depend on the factorization of P that is employed 
(Vidyasagar,1985, p. 1 10) [14]. 

Let h ( p,q ) be the impulse response sequence that corresponds to H ( P,Q ). Let also 
r, = { '«(*) ) *=o . * = 1 > 2, 3 and q = { q(k ) } " =0 , where ti(k) and q(k ) are real matrices of 
appropriate dimensions. Then, the sequence h (p,q)= {h(k) }*=o is given by : 


h(k) = 


h lid) h u {k) 


’ ' hmn ) 


t,{k)~ £ 

j=° 


t 2 (k-j) £ q(X)t 2 (J-^) 

X=0 


( 1 ) 


Similarly, let s i ,s 2 ,s 3 ,s 4 be the impulse response sequences, members of M (/ 1 ), that 
correspond to Sj ,S 2 ,S 3 ,S 4 respectively and s,= { S[(k) }” = o i= 1,2,3, 4. Then, the quadratic 
constraint is satisfied iff all the elements of the impulse response sequence that coresponds to the 
LHS are equal to zero. Thus the following infinite set of quadratic equality constraints is 
obtained: 


k 

I 

j = o 


\ q(J)s 2 (k~j) + s 3 (k-j)q(j) + q(k-j ) 

is 4 (j-^)qU) 

>1 

h=0 JJ 


= 0,V*>0 (2) 


The objective function of (DPP) becomes: 


\\H(P,Q)\\^ = \\h(p,q)\\r= max £ £ 

i=l. j = l k=0 


ti(k)-Z 
y=0 


t 2 (k-j) X q(k)t 3 {J-X) 

x=o 




where 


kj(k)=\ r,(*)- I 

7=0 


tiik-j) ^qO^Hij-X) 

x=o 




Based on (1) and (2) (DPP) is transformed into the following constrained / 1 -optimization 
problem: 


subject to. 


v = inf |/i (p,q)\\i l 


fk(q) = 0 , V*> 0 


(DPPs) 


For a given sequence q and a value of k,f k (q) is an mxn matrix with entries: ft J (q ) , , 

j=l,...,n. Then (DPPs) can be reformulated as: 
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(DPPs-a) 


v = inf \\h(p,q)\\r 

subject to, 

m n oo 

lit >fi’(.Q)l=0 

i= 1 j- 1 *=0 

These optimization problem are infinite dimensional; the optimization variable (q) lies in 
an infinite dimensional linear space and the constraint is also infinite dimensional. In the follow- 
ing, it will be shown that one can obtain solutions, arbitrarily close to the solution of these prob- 
lems, by solving appropriately constructed finite dimensional optimization problems. 

To reduce (DPPs) to a finite dimensional optimization problem two types of truncations are 
performed: 

Truncate the constraint 

V w = inf \\h(p , (?)||/i (DTI) 

subject to, 

fk(q) = 0 , k= 


and 


Truncate the variable q and the constraint 

Vn = inf \\h(p,q)\\i' (DT2) 

subject to, 

fk (q )=o , k =0, - • • ,M 

where c 4>o is the set of all l l mxn sequences with their first N + 1 elements nonzero. 

In the following, it is assumed that (DPPs) is feasible, thus feasibility of (DTI) is 
guaranteed. Under this condition, feasibility of (DT2) is guaranteed provided that N >M. The 
following theorem establishes the relationship between (DTI) and (DT2). 


Theorem 1 


lim vjy = v w 

jV— 


Before proceeding with the proof of the theorem, let us consider the following sets: 
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Q <= Imxn 


I fi(q)= 0 , i = 0,---,M 


^Al- 


and i 

? 

G 0 m = j<7i 6 l/;(<7i)=0, i=0,---,M>cG m . j 

The first set is the feasible set of (DTI). The feasible set of (DT2) is a subset of the second set. 

Then the following lemma establishes a result that will be used in the proof of theorem 1. 

i 

I 

Lemma 1 

The set G 0< m is dense in Gm- , 

Proof of Lemma 1 1 

To prove that the first set is dense in the second it suffices to show that : given 5 > 0 and I 

q e Gm there exists q \ e G q m such that \\q - q \ ||/> < 5. 

The set of M + 1 constraints f l (q) = 0 , i= 0, ■ - • ,M involves only the first M+1 elements of 
the sequence q. Any other sequence with the same M+1 first elements satisfies the constraints. 

Therefore, q can be partitioned as : q = f q a I ] e Gm, where q a e R w+t is the vector of ele- 1 

ments of q that appear in the constraints. Clearly, qt, e / 1 and there exists a sequence q \ ^ e <|> 0 j 

such that for given 5 > 0 => \\q b - q i <b ||/> <8. The sequence q\-[q a I q 1 is a member of 
G 0 'M ar >d : • 

Ik — <7 J II/ 1 =ha-Qah +hb-qi,b\\l' < 5 - 

□ i 

Now, the proof of theorem 1 follows. 

Proof of Theorem 1 1 

We want to prove that lim = v w . 

n -*° o : 

Equivalently, for a given e > 0 we want to prove that there exists N such that Iv** - |< £. ; 

Since, v M = inf || h (p,q) ||/i , for a given E > 0 there exists a q e G M such that : : 

v"-|*<M>li'|<f ( 2 ) ; 

Since \\h (p,-)!/ 1 is continuous on Gm then, given qe Gm there exists 8 > 0 such that for all 
q'<=G M , 

\\q-q'\\i' <5 => \\h(p,q)\\i' -\\h(p,q')\\^ <-| 

Since G 0i m is dense in G M => there exists q \ e G q< m : ||<7 -<7 i II/ 1 <8. Let N be the smallest 
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number such that q \ e $ 0tN and q \ satisfies the last relationship. 
Consequently, 


\\h(p,q)\\t' -\\h(p,q\)\\i 

From (2) and (3) the following relation is obtained : 


E 

< 2 


( 3 ) 


< E 


Since < vjv and 0<vjv < ||/i (p, <7 1 ) ||/ 1 we finally obtain : v^ — v^f|<£. 


□ 


This theorem establishes the first approximation result. However, its application assumes 
that for any given N and M the value V/v is known, or can be computed. The solution of (DT2) 
involves the minimization of an infinite sum. 

Given an l l mxn sequence £, its t ] (L) sum is defined as : 


\m i'(L) 


max £|&J/‘(L) 

i=i. • ■ • ,m J= j 


n 

max £ 
i=i. 


X IM*)I 


*=0 


Define the set of sequences q e <[) 0i /v that are norm bounded by some positive number B : 

q e <> = | q e 4> 0 ,a/ ; Ik II / 1 ^ B 

The set <|> J N is finite dimensional (so is 0 OiA /). closed, bounded and therefore by the Heine-Borel 
theorem it is compact. Consider a formulation of (DT2) where the variable is norm bounded. 

vjv(B) = inf \\h(p,q)\\i' (DT2a) 

<?6 0®iV 

subject to, 

fi(q) = 0 , i =0, • • • ,M 

Clearly, vjv ^ vjv (B) for all values of B. 

Let q e be a suboptimal solution of (DT2), i.e.: 

V e > 0 , 3qe <|) 0 jv ; 0 < \\h (p,qji\i' - vjv < £ 

The norm of q is finite. Let |k|k = B (e). Existence of the solution of (DT2a) is guaranteed by 
the compactness of its feasible set (see proof of lemma 2). Let q be the solution of (DT2a). 
Then, q' satisfies the relation: 

vjf — Vn (B (e)) — 1| /i ip,q') I/ 1 <\\h{p,q)\\o 

Then, combining these statements the following is obtained: 
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V e>0, 3B (e) ; 0<vjv(5 (e))-vjv <e 


If B is selected to be sufficiently large, the solution of (DT2a) is arbitrarily close to the solution 
of (DT2). In subsequent sections it will be demonstrated that the calculation of bounds for the 
optimization variables is feasible. 

Based on this discussion, the following optimization problem is formulated: 

Vn,l - inf ||/t(p,<7)||/ l (L) ( DT3 ) 

subject to, 

/*(<?)= 0 , k = 0 , ■ ■ ■ ,M 

The following lemma establishes the fact that the solution to (DT2) can be obtained through 
iterative solution of (DT3) for increasing values of L. 

Lemma 2 

lim v# L = v#(S) 

L — 

Proof 

The sequence { vj^ } 2T=o is nondecreasing and bounded: 

VL> 0 : 0 <vft L <vft(B). 

As a result : 3a > 0 ; lim v)vx = a. 

This limit cannot be greater than vj v (B). Assume that a < vj tf(B). 

The feasible set of (DT3) is compact. Indeed, it can be written as the intersection of a finite 
dimensional, closed and bounded set (§o,n) with a finite dimensional closed set: 

Gq M — §o,N <7 1 S tyo.N i fi ( Q 1 ) ~ ® f — 

As a result the feasible set is finite dimensional, closed and bounded. 

From the continuity of the / 1 (L) norm, and the compactness of the feasible set it follows 
that (Luenberger, 1969; p. 14) [ 8 ]: 

Vl>0,3<7 L g r\G oM ;v%' L =\\h(p,q L )\\i' (L) 

Compactness of the feasible set implies that the sequence { qi 1^=0 has a subsequence 
{ Q Lk } r =0 that converges in <$>% O C 0 , M : 

lim <74 = <7 e <t >o,N O G 0 'M 

Then, 
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limv^f^= lim \\h (p,<7z.)||/ 1 (L) - bm WhiP^L^Wi'm ) _ a 

l^oo * L-»~ 


This contradicts the assumption that a < v^f (B). Hence, a - v# (B). 


□ 


In summary, lemma 2 establishes that solution of (DT3) for increasing values of L identifies the 
solution of (DT2a). Based on the "equivalence" of (DT2a) and (DT2), it can be said that this 
iterative procedure identifies e - optimal solutions of (DT2). 

Remark : In the case where the coprime factors of P have been constructed to be FIR’s, for each 
q G <]) 0 N the sequence h ( p,q ) will have a finite number of nonzero elements. This number is 
known and depends only on N. Therefore the norm of h ( p,q ) can be exactly calculated by a 
finite sum. In this case, the exact solution of (DT2) is obtained in one step. 

The relation between the two different types of truncated problems has been established 
through Theorem 1. In the remaining of this section the connection between (DTI) and (DPPs) is 
shown. For a non-negative real number 5 consider the following sets: 


m n 


G( d)=\qei 


mxn 


I 

;=i j 


I I 

1 k = 0 


fi J (q) I <8U M (5)H<?el, 


mxn 


M 

m 

i=l j=\ k= o 


A j (q) i 


Then, the following theorem holds: 


Theorem 2 

Assume that G (•) and G M ( ) represent upper semicontinuous mappings from the non- 
negatives reals to subsets of / 1 . Then, 

lim v M = v. 

M — 


Before we proceed with the proof of Theorem 2, two intermediate results will be presented. 
Consider the two optimization problems: 

f m n °° 

v(co) = inf i II^J (P><7)||/' + ^ ^ 

i=] J =1 k =° 


{ m n M 

11^ (p>q)\\i' X 2 

i=l ;=1 *=0 J 

It is easily verified that (PI) is the penalty function formulation of (DPPs-a) (Luenberger, 1969, 
pp.302-305) [8]. For this type of problems the following lemma can be shown to hold: 

Lemma 3 

Let {o/}r=o be an icreasing sequence such that lim o/ = °°. Let also G (•) be an upper semicon- 

r — 

tinuous mapping from the non-negative reals to subsets of /i, X n- Then, the following statements 
hold: 


I fk } (q) I 


fH J (q) I 


(PI) 


(P2) 
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1. v(co r+1 ) > v(co r ) 

2. v>v(co r ) 

3. limv(co r ) = v. 

r — >oo 

Proof 
Part 1 & 2 

The proof of Part 1 & 2 is performed in a similar way as the proof of statements 1 & 2 in Lemma 
1 in Luenberger (1969, p.305) [8]. 

Part 3 

From Part 1 & 2 it follows that the limit exists and is less than or equal to v. The proof is similar 
to the proof of Part 3 in Lemma 4.2 in Sourlas and Manousiouthakis (1992) [13]. 

□ Employing 

the same line of arguments the following corollary can be shown to hold: 

Corollary 1 

Let G m (- ) be an upper semicontinuous mapping from the non-negative reals to subsets of Imxn- 
Then, 

1. V w (0)'' +1 ) > v M (to r ) 

2. \ M > v M (ci/) 

1. limv A/ (o)' r ) = v M . 

r— » oo 

For a given value of to, Lemma 4 establishes the relation between v^(co) and v(co). 

Lemma 4 
lim v M {(£>) = v(to). 

M ->« 


Proof 

The proof is based on the fact that the sequence {/*(<?) } T=o belongs to /^,xn thus allowing the 
approximation of the infinite sum in the penalty term in (PI) with a finite sum. The detailed 
proof can be found is similar to the proof of Lemma 4.3 in Sourlas and Manousiouthakis (1992) 
[13], □ 

Now the proof of Theorem 2 follows. 

Proof of Theorem 2 

The proof consists of two parts. First existence of the limit will be shown and then conver- 
gence to v will be demonstrated. The sequence v M is nondecreasing, and bounded above by v. 
As a result, it converges to a real number: 

~ hi 

a= lim v < v. 

M —4 oo 
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( 4 ) 


Using Corollary 1, Lemma 3 and 4 the last relation can be rewritten as: 

a= lim v w = lim lim v M (oo)<v= lim lim v^co) 

A/— >oo Af — >oo (o— >oo ( 0 — >°° M >°° 

In the remaining of the proof it will be shown that the two iterated limits are equal. 
Using the proof technique of Lemma 3 it can be established that: 

v M (o))< limv M (co) , ViW.Vco 

CO — >°° 

This implies that: 

v(oo)= limv M ((o)< lim limv M ((o) , V co => 

M — >°° M — >o° co— 


v= lim limv^fto)^ lim limv w (0)) = a 

(0—^03 /V/^OO M — > OO (0 — >00 

In view of (4), the last relation implies that a= v. 

As a result, the following statement has been proven: 


a = lim Vm = lim lim Vyv^ = v 

M — > oo M — ) oo /V->~ 


□ 


4.2. Computational Procedure 

In view of theorems 1,2 and lemma 2, the value v is obtained through a limiting procedure, 
which involves solution of (DT3) for increasing values of L, N and M. The finite dimensional 
optimization problem can be formulated as a nonlinear program, which can be solved by finite 
dimensional optimization techniques. For a particular example, the structure of the nonlinear 
program is given in appendix A. 

The computationally intensive part of the procedure is the solution of (DT3), which is a 
nonlinear programming problem, for different values of L, N and M. The resulting nonlinear 
programs are nonconvex, due to the existence of the quadratic equality constraints. Nevertheless, 
global solution of these problems determines the globally optimum performance that can be 
achieved by a certain decentralized structure. If (DT3) is solved locally, then the limit of the 
resulting sequence will identify an upper bound to the / 1 optimal decentralized performance. 

4.3. Global Optimization Approach 

Global optimization approaches that solve the general nonconvex optimization problem 
have recently been developed. Manousiouthakis & Sourlas (1992) [10], presented a procedure 
that is based on the transformation of the original nonlinear optimization problem into one that 
has convex constraints and objective with an additional separable, quadratic, reverse convex 
constraint. Employing this transformation procedure, (DT3) becomes a convex programming 
problem with an additional reverse convex, quadratic and separable constraint. This problem 
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can then be solved and its global optimum can be identified through the use of a branch and 
bound type of algorithm. 

The implementation of this global optimization algorithm requires the existence of valid 
upper and lower bounds for all elements of the the sequence q that appear in the quadratic equal- 
ity constraints in (DT3). Bounds on these variables can always be obtained from local minima 
information, namely the value a local minimum of (DT2) for fixed values of N and M. One can 
always obtain such information if the solution of (DT3) for fixed N, M and increasing L is per- 
formed using local optimization algorithms. The generated sequence of values converges to the 
value of (DT2) at a local minimum. Let 8 i ocal be this value. This is an upper bound to the glo- 
bally optimal value of (DT2). Thus, for q e and for each value of L, it holds that: 

n L 

\\h(p,<l)\\t'(L) * ?>iocai <=> max £ Z \hij(./c)\<8 loca i => 

1=1, • • • ,my = 1 k=0 


hij(k) 


&local * S 'local 


,/=!,•*• , k=0, - ,L 


(5) 


Then, one can obtain lower (upper) bounds on all elements of the sequence q through the solu- 
tion of a minimization (maximization) problem with objective the corresponding element of the 
sequence and constraints the inclusion realtions that appear in (5). These optimization problems 
are linear. 


5. Example 

In this section the computational procedure introduced in section 4.2 is applied to the fol- 
lowing 2x2 example. 


P(z) = 


3 z + 1 

z 

z 2 + 2.5z + 1 


0.5 z +2 
z 

3z + 4,5 
z 


The objective of our analysis is to determine the best possible achievable performance for 
a decentralized control system featuring the pairings ( (yi,«i );(>’2> M 2) )■ 

Parametrization of all Decentralized Stabilizing Controllers 

The process P is stable. Therefore one can select coprime factors as follows: 

Np = Np = P , D P = Y = I , D P = Y = I , X=X = 0 


According to (2), sec.3.1, any stabilizing controller for P is parametrized as: 


c=q(i-pq r 1 


011 012 
021 022 


€ M(S) 


( 1 ) 


The controller C is decentralized iff: 

0 21 ^ 12 =P 210 12 
0 12 =“^ 12(0 11 022 “0 12021 ) 
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The Input - Output map between disturbances and output is given by: 

H(P,C):u 2 ->y 2 -,H(P,C)=(r + PCr'P = U-PQ)P 
Problem Formulation 

The disturbance rejection problem is formulated according to the guidelines introduced in 
section 4. 1 . First, the known bounds on the disturbance and the desired bounds on the objective 
are defined: 

Disturbance: «2(0 e [-1, 1] x [-1, 1] , Vr 

Output: y\(t)& [-1,1] x[-l,l] , Vr 

Then, (DPPs) is readily formulated and transformed to (DT3) according to the procedure 
presented in Section 4. Based on Appendix A, (DT3) is in turn transformed into a nonlinear pro- 
gramming problem. Since the coprime factors are F.I.R. then the solution to (DT2) is obtained 
in one step. For N=M=0, 1 ,2 the globally optimum value of the corresponding optimization 
problem has been identified. The complete sequence of values that converges to the value of v 
for this particular problem is shown in the following table: 


N 

M 

— z ® — 

v/v 

N 

M 

— 

vn 

0 

0 

8.00 

5 

5 

5.20 

1 

1 

6.70 

6 

6 

5.16 

2 

2 

5.85 

7 

7 

5.13 

3 

3 

5.51 

8 

8 

5.13 

4 

4 

5.32 

9 

9 

5.13 


Hence, with e= 10 -2 the value of the / 1 - optimal decentralized performance problem was found 
to be 5.13. When only the linear constraints are considered the value of the global lower bound 
to the Z 1 decentralized performance is 4.98. The /' - optimal centralized performance, for the 
same set of specifications was to found to be equal to 4.72. 


6. Conclusions 

Based on the Manousiouthakis parametrization of all decentralized stabilizing controllers, 
the Z 1 optimal decentralized performance problem has been formulated as an infinite dimen- 
sional optimization problem. This problem was transformed into a finite dimensional one 
through the introduction of appropriate truncations. Theorems that establish the equivalence (in 
the limit) of the original problems to their finite dimensional approximations were proven, and a 
computational procedure was proposed. It has been established that solution to the optimal 
decentralized performance problem amounts to global solution of a series of quadratically con- 
strained programming problems. If locally optimal solutions are identified for each of the finite 
dimensional problems, the limit of the corresponding sequence of values will be an upper bound 
to the optimal Z 1 decentralized performance. Based on this work, one can actually evaluate the 
best performance achievable by a given decentralized structure. 
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Appendix A 

The optimization problem (DT3) can be formulated as a nonlinear programming problem. 
The steps that make this transformation feasible follow. For illustration purposes the 2x2 case, 
with stable plant, the same as the one in the example, is considered. From the definition of the 
/’(L) sum : 
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L 

h(p,q)e 1 2x2 => \\h(p,q)\\i'(L) = 5= max £ |Aii(*)|+l*/2(*)| 

,= 1'2 1=0 


Using the definition of the maximum as the least upper bound, (DT3) is finally transformed 

into: 


v nl = inf 5 

.,7=1.2 

k=0, ■ ,L 


(A.l) 


subject to, 

k 

I 

i=0 


/*(?)= X [?12(OP2l(*-O-<72l0>12(*-O =0 
; L 


q 1 1 W)q 22(1 ~j) — Q n(j)Q 21 O'— y) 


£*(<7)=<7i2(*)+ X/U2(*-0X 
.=0 y=0 L 


q,j(k) = Q , /, j = 1,2 , A/+1 < /: <L 
-5 iy (it)</iy(it)<8y(it) , i, j= \ , 2 , k = Q, • ■ • ,L 


k = 0, • ■ ■ ,M 


= 0 


L v 

X 5,1 
*=0 L 


(k) + b l2 (k) 


<5 , i'=l,2 


where fidq) and gk(q) are scalar constraints resulting from the application of the parametrization 
to this particular case (2x2 controller, stable plants). 
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